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Motivation 


Introduction and Motivation 
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Motivation 


Path integral with complex weights appear in many important 
physics: 


e Finite-density lattice QCD, 
spin-imbalanced nonrelativistic fermions 


ə Gauge theories with topological 6 terms 


ə Real-time quantum mechanics 


Oscillatory nature hides many important properties of partition 
functions. 
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Example: Airy integral 
Let's consider a one-dimensional oscillatory integration: 


a= | i T as 
p 3 - 


RHS is well defined only if Ima = 0, though Ai(z) is holomorphic. 


Horo nap 


exp( osx ested) 


AA AAAI 


-8 6 4 -2 0 2 4 


Figure: Real parts of integrands for a = 1 (x10) & a = 0.5i 
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Contents 


ə How can we circumvent such oscillatory integrations? 
=> Lefschetz-thimble integrations 
[Witten, arXiv:1001.2933, 1009.6032] 


ə Applications of this new technique for path integrals 
> Quantum tunneling in the real-time path integral 
[YT, Koike, Ann. Physics 351 (2014) 250] 
> Phase transition of the matrix model 
[YT, arXiv:1412.1891; Kanazawa, YT, arXiv:1412.2802] 


Yuya Tanizaki (University of Tokyo, RIKEN) Lefschetz-thimble path integral Jan 9, 2015 @ OIQP 


5/2 


Lefschetz-thimble integrations 


Introduction to Lefschetz-thimble integrations 
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Again Airy integral 


Airy integral: 


3 
Ai(a) = f T epi (5 + ar) 


The integrand is holomorphic w.r.t x 
=> The contour can be deformed continuously without changing the 
value of the integration! 
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Lefschetz-thimble integrations 


Steepest descent contours 
What is the most appropriate contour for our purpose? 


Re[iS(a,a)] should be made as small as possible. 
< The contour should be perpendicular to Re[iS(, a)] = cosnt. 
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Lefschetz-thimble integrations 


Steepest descent contours 
What is the most appropriate contour for our purpose? 


Re[iS(a,a)] should be made as small as possible. 
<> The contour should be perpendicular to Re[iS(, a)] = cosnt. 


Steepest descent ones J must satisfy Im[iS(x,a)] = const. 
because of the holomorphy. 


m 


if 


xp +i0.1. 
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Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (o = 1,2) for the Airy 


integral: 
2: of Zepi (Z +a). 


No: intersection number of the steepest ascent contour Ke and R. 
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Figure: ee ibis J and duals K (a= = 'exp(0. 1i), exp(mi)) 
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Generalization to multiple integrals 


Model integral: 


z= | dx, ---dz, exp Z(x;). 
Rr 
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Lefschetz-thimble integrations 


Generalization to multiple integrals 


Model integral: 


z= | dr): - d£n exp Z(x;). 


What properties are required for Lefschetz thimbles J? 
@ J should be a n-dimensional object in C”. 
© Im[Z] should be constant on J. 
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Lefschetz-thimble integrations 


Short note on technical aspects 
Complexified variables (a = 1,..., n): Za = Za + ipa. 
Regard x, as coordinates and p, as momenta, so that Poisson 
bracket is given by 


N [seat - so 


a=1,2 


O1 Opa OLa Opal 
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Lefschetz-thimble integrations 


Short note on technical aspects 
Complexified variables (a = 1,..., n): Za = £a + ipa. 
Regard x, as coordinates and p, as momenta, so that Poisson 
bracket is given by 


u Of Og ög Of 
gh = 2 Ee Ops: On, an 


Hamilton equation with the Hamiltonian H = Im[Z(z,)]: 


df (x, p) dza OT 
ig eee (= dt --(#)) 


This is Morse’s flow equation: “Re[Z] < 0. 


= We can find n good directions for J around saddle points! 
[Witten, 2010] 
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Multiple integral: Lefschetz-thimble method 


Oscillatory integrals with many variables can be evaluated using the 
“steepest descent” cycles Js: 


d?r eiS(2) = ne | d"z ee, 
I 2 Jy 


Js are called Lefschetz thimbles, and Im[iS] is constant on it. 


Ne. intersection numbers of duals X, and R”. 
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Searle E EEC Harmonic oscillator 


Example of QM: Harmonic oscillator 


The action functional 


Saddle-point condition (=Euler-Lagrange eq.): 


d?z.ı(t 
zal ) = — zalt). 


The classical solution: 
Te — T; cos T 


zalt) = ao ae sin t + x; cost. 


Lefschetz-thimble path integral 
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bent caeli 
Example of QM: Flow equation 


Flow equation 


O aala) A ee 
—Az(t;u) = —i | — z(t; u), 
Ou ot? 

Boundary conditions: Az(t; —co) = 0, and Az(0; u) = Az(T, u) = 0. 
The set of solutions Az are spanned by 

eir/4 exp (e — 1) u sin”, (nn>T), 
e-ir/4 exp (i — (7) u sin”, (nt <T). 


Azn(t; u) = 
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Penit eee 
Example of QM: Lefschetz thimble 


Lefschetz thimble J (~ R”): 


. net = wet 
—im/4 Bee im/4 MR 
za(t) + € X usin T +e ` arsin = | ae € R 
é=1 f=v4+1 
v: the maximal non-negative integer smaller than T/r. 
Integration measure on J becomes 
V . (oe) . (oe) 
/ Dz = nf I] e idan I] ei dan = eeN |T] Vidap. 
I n=1 m=v+1 f=1 
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binet sa. 
Example of QM: Path integral 


Feynman kernel for the harmonic oscillator: 


| 
© 


TZ 
E Zei] _ N vin E 
Tze) _ ‚TU 1 1 
xp (A E / Var jl T/r] 
= Alea oe us 
u 2rihl aA 2 


| 
© 


Appearance of the Maslov-Morse index v is quite manifest. 
[YT, Koike, Ann. Physics 351 (2014) 250] 
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Haven oral es 
Proposed Applications of Lefschetz-thimble methods 


ə Analytic properties of partition function for CS theories 
[Witten] 


ə Real-time description of tunneling 
[YT, Koike; Cherman, Unsal] 


oe Sign problem of the lattice Monte Carlo simulation 


[M. Cristoforetti, et al., Kikukawa, et al.] 


e@ Geometrization of resurgence transseries 


[Basar, Dunne, Unsal; Cherman, Dorigoni, Unsal] 
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Quantum tunneling 


Quantum tunneling and real-time path integral 
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Quantum tunneling 


Double-well potential 


dz 2 2 2 
= — 12 = 
=) + (z°—1)"=p 


Two different origins of oscillations 
=> Solutions show double-periodicity! 


Complex-energy conservation: ( 


real-time oscillation imaginary-time oscillation 
A 


>x 


The list of parameters p can be obtained as (n,m € Z) 


pd PT) aa NOS E (ste: 


v2p v2p 2 
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Quantum tunneling 


Complex classical solutions 
Analytic complex solution with z(0) = 2;, z(T) = xy: 


p2—1 = 2p 1l+p l+p 
= 2 i; ’ 5 
2(t) = | 3p sd (v pt + sd (= 3p DD 
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Quantum tunneling 


Complex classical solutions 
Analytic complex solution with z(0) = x;, z(T) = xy: 


. jal 1 2p lt» 1+p 
z(t) = 3p sa ( V/2pt + sd (= op) oe 


List of possible parameters p: (YT, Koike, Ann. Physics 351 (2014) 250) 


Imp) 
150F 


. ` ook 7 


Re(p) 
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Quantum tunneling 


Complex classical solutions 
Analytic complex solution with z(0) = x;, z(T) = xy: 


=ï 7 2p I+p I+p 
dl ./2pt + sd! j ae, 
2p ( a (eta VE) 2p 


List of possible parameters p: (YT, Koike, Ann. Physics 351 (2014) 250) 


p? 


Im(p) 
150F 


. ` ook 7 


Re(p) 


=100F 


= 150 


Let’s describe quantum tunneling semi-classically using these 
“complex” solutions. 
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Quantum tunneling 


Imaginary-time instantons “=” Sphalerons 
Instantons in imaginary times (p œ 0): 


(n,m)=(0,0) (n.m)=(1,0) 
\ Wick rotation 
Sphalerons in real times (p ~ 1): 
(n,m)=(0,0) = (n,m)=(1,0) 
| Fa 9 
p am 4 EJ 8 w 
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Quantum tunneling 


Important! 


Any complex solutions behave as sphalerons in the limit 
Zee. 


Energy p? of the solution is determined by 


nos (p+ 1)/2p) | iK(y (p — 1)/2p) L pis tere: 


m = 
y 2p y 2p 2 
Labels n,m roughly determines the number of oscillations during T. 
=> If n,m < T, the solution must lie on the top. Thus, p ~ 1. 


(n,m) = (2,1) 
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Quantum tunneling 


Complex solutions for quantum tunneling 


In order for the one-instanton action iS[z] ~ —Sinst. = —4/3, 
complex solutions must be (highly-)oscillatory in the complexified 
space, i.e., n,m ~ O(t; — ti) 

(YT, Koike, arXiv:1406.2386; Cherman, Ünsal, arXiv:1408.0012): 


(n,m) = (31,30) 
15 (n,m) = (52,50) 


AA a) 


w 


i 80 
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Summary & Perspectives 


Summary & Perspectives of this section 


Summary 


ə Real-time dynamics becomes calculable in a nonperturbative way 


with Lefschetz-thimble path integrals. 


ə Exact semi-classical description of quantum tunneling is 
considered; 
@ Imaginary-time instantons are connected to sphalerons by 
analytic conti. 
@ Real-time tunneling must be described by highly-oscillatory 
complex solutions. 


Perspectives 


ə Application to real-time dynamics of QCD: Sphaleron rate r. 
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Gross-Neveu-like model 


Phase transition associated with symmetry 
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Gross-Neveu-like model 


Motivation 
What is SSB in this context? 


Let’s consider two questions using zero-dimensional fermionic models 
with N species: 


ə What happens when discrete symmetry is broken at large N? 


ə Can we also treat continuous symmetry? 
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Gross-Neveu-like model 


0-dim. Gross-Neveu-like model 


The partition function of our model study is the following: 


Zu(G,m) = | apaw exp [Sua + mya t ur ( ya) | l 


The Hubbard-Stratonovich transformation gives 


N 
Zn(G,m)=4/ “a faces, 


— log[p* + (o + m)’]. 


with 


For simplicity, we put m = 0 in the following. 
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Gross-Neveu-like model 


Properties of S(c) 
S has three saddle points: 


_ OS(z) _ 22 2z 


0 u 
Oz G P+2 


= z=0 +/G-P. 


Figures for G = 0.7e*°" and p? = 1 [Kanazawa, YT, arXiv:1412.2802]: 
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Figures for G = 1.1eF®"! and p? = 1 [Kanazawa, YT, arXiv:1412.2802]: 
8 P 
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From these figures, we learn that, for real G, 


ə z = is the unique critical point contributing to Z if G < p°. 
ə All three critical points contribute to Z if G > p°. 
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Gross-Neveu-like model 


Stokes phenomenon 


The difference of the way of contribution can be described by 
Stokes phenomenon. 

< At some special values of coupling, several critical points are 
connected by the flow. [Witten, arXiv:1001.2933, 1009.6032] 


Figures of G-plane for ImS(0) = ImS (z4) [Kanazawa, YT, arXiv:1412.2802]: 


-1.5 


“10-05 00 05 10 15 20 


RG 
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Gross-Neveu-like model 


Dominance of contribution 


The Stokes phenomenon tells us the number of Lefschetz thimbles 
contributing to Zy. 


However, it does not tell which thimbles give main contribution. 


In order to obtain (a) # 0 in the large-N limit, z+ should dominate 
z=0. 


> ReS(z;) < ReS(0) 
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Connection with Lee-Yang zero 


Blue line: ImS(z+) = ImS(0). 
Green line: ReS(z+) = ReS(0). 
Red points: Lee-Yang zeros at N = 40. [Kanazawa, YT, arXiv:1412.2802] 
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Gross-Neveu-like model 


Conclusions for studies with GN-like models 


@ Decomposition of the integration path in terms of Lefschetz 
thimbles is useful to visualize different phases. 


@ The possible link between thimble decomposition and Lee-Yang 
zeros is indicated. 
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Lefschetz thimbles with continuous symmetry 


Lefschetz thimbles with continuous symmetry 
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uri erensia 
Lefschetz thimble with approximate O(2) symmetry 


Let's consider approximately O(2)-symmetric sigma model: 
S(@) = -(a? — 1)? +eoı. 


“Angular momentum” gives an extra conserved quantity at € = 0! 


Po = 1201 — 102. 
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Lefschetz thimbles with continuous symmetry Simple bosonic model 


Lefschetz thimble with approximate O(2) symmetry 


Let's consider approximately O(2)-symmetric sigma model: 
Slo) = —(@” — 1)? + e€. 
“Angular momentum” gives an extra conserved quantity at £ = 0! 


Po = M201 — 102. 


At £ = 0, the flow is singular, and naive construction of J fails 


down. 
Under this circumstance, how does the flow look like? 
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Lefschetz thimbles with continuous symmetry Simple bosonic model 


Structure of the flow equation 


Hamilton equation (= flow equation): 


F 2 
al) +ecos6, 
E 


dt 
SLR 1 — 3r? + Dp — . - sin, 
dt 

dd 2p, Po sin 9 

— = € ’ 

dt r T 
dpo _ : Po 
Fr (p, sind + 2 cos 8) ; 


In the limit € — 0, flows along symmetric direction disappears. 
=> Naive construction of Lefschetz thimbles breaks down at £ = 0. 
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Lefschetz thimbles with continuous symmetry Simple bosonic model 


Flow with approximate O(2) symmetry 


Set of saddle points ate = 0: o? - n? =1l&o-n=0. 


Flows in the oj09-plane & in the o17-plane for small e: 


ə Due to the approximate O(2) symmetry, flow walks around the 
pseudo saddle points. [YT, arXiv:1412.1891] 
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Lefschetz thimbles with continuous symmetry 


Flow with approximate O(2) symmetry 
Schematic pictures for Lefschetz thimbles of the system 


(o° -— 1) + en) 


po/r oy =rsind 
A 


7 


J1-2/2.0) 


01 =rcosd 


Sare/2,n) 


=> In order to get an O(2) symmetric result, we must sum up two 
thimbles Ji+e/2,n) and Ja-e/2,0)- [YT arXiv:1412.1891] _ 
Yuya Tanizaki (University of Tokyo, RIKEN) 


Lefschetz-thimble path integral 


DA 


Lefschetz thimbles with continuous symmetry Simple bosonic model 


Singularity of symmetry-restoration limit 


What is the fate of Z14</2,n) and Ja-—e/2,0) in the limit € — 0? 


Integration on Jı1-./2,0): 


/ d’zexp [—S/h] 
JTAa-e/2,0) 


~ -i L AdA exp CA m dd exp (-2 cosh 6) 
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Lefschetz thimbles with continuous symmetry Simple bosonic model 


Singularity of symmetry-restoration limit 


What is the fate of Z14</2,n) and Ja-—e/2,0) in the limit € — 0? 


Integration on Jı1-./2,0): 


/ d’zexp [—S/h] 
JTAa-e/2,0) 


~ -i m AdA exp (5) dd exp = cosh 6) 


Logarithmic divergence: f dọ exp (- € cosh ¢) ~ 21n £ 


= Thimbles related by the symmetry need to be a up to give 
finite results! 
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Lefschetz thimbles with continuous symmetry Nambu—Jona-Lasinio-like model 


Nambu-Jona-Lasinio-like model 
Let's apply this observation to the following NJL-like model: 


N 
Zua) = [asa exp (>: Wa(ip +m) 


N ~~ N 2 
(> int) ~- p inh] 
a=1 a=1 


Bosonization: 


Zua) = oe exp (—NS(o,7)) , 
TG 
with 
2 2 
Oo +7 
S(o,m) = — log (p + (o +m)? +977) + G 
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Lefschetz thimbles with continuous symmetry 


Lefschetz thimbles for NJL-like model 
In this case, m(# 0) breaks the O(2) symmetry: 


+ or 
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Lefschetz thimbles with continuous symmetry  Nambu—Jona-Lasinio-like model 


Conclusions for the continuous symmetry 


e Lefschetz thimbles with approximate continuous symmetries are 
computed. 

ə The remnant of the symmetry implies existence of slow flows 
on the set of pseudo saddle points. 

e Such slow flow relates different thimbles via the approximate 
symmetry. 
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